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We study the dynamical properties of recently introduced 
frustrated lattice gas models (IFLG and Tetris) for granular 
media under gentle shaking. We consider both the case where 
grains have inter-grain surface interactions and the case where 
they have not, corresponding, for instance, to the presence or 
absence of moisture in the packs. To characterise the grains 
packing structure, we discuss the properties of density distri- 
bution. In particular, we consider the phenomenon of grains 
domains formation under compaction. New results amenable 
of experimental check are discussed along with some impor- 
tant differences between the dynamics of the present models. 



I. INTRODUCTION 

Non-thermal disordered systems as granular media 
show a variety of dynamical phenomena with interesting 
properties For instance, important physical fea- 

tures were recently discovered in shaking experiments, 
where the dynamics of grains is driven, under the effects 
of gravity, by sequences of shakes ( "taps" ) ^ . A very 
common phenomenon recorded in those experiments is 
compaction: when a box filled with loose packed sand is 
shaken at low amplitudes, T^x {^ex — a/ff, where a is the 
shake peak acceleration and g the gravity acceleration 
constant), the bulk density of the system increases with 
the number of shakes. Interestingly, the Chicago group 
has experimentally shown that compaction is logarith- 
mically slow H]. Several models have been proposed to 
describe these kind of dynamical behaviours P-p^ , but 
a detailed investigations of its microscopic nature is still 
incomplete. 

In the present paper we discuss some aspects of such a 
problem by focusing on the details of the packing struc- 
ture. Our analysis is developed in the context of re- 
cently introduced lattice gas models for granular me- 
dia the IFLG and Tetris. These models are 
very simple and thus very schematic, but they interest- 
ingly describe, in a single framework, many dynamical 
properties of granular media ranging from logarithmic 
compaction "irreversible-reversible" cycles and aging 
PI , anomalous dynamical responses [p^ , to segregation, 
avalanches effects and several others ]13[ | . As a general- 
isation of the q-Model they have been also studied 
to describe "scalar force" patterns. 

To understand the grain packing structure, we 
recorded the bulk density distribution in the present 
models, which in some cases turns out to be definitely 
non Gaussian. We propose a variant of our models to 



be able to describe also grain-grain surface interactions, 
which are to be considered for instance in cases of charged 
grains or non dry media. Interestingly, we find that the 
dynamics of compaction is not qualitatively affected by 
these details of grains interactions, at least up to when 
the coupling strength is of the order of the effective tem- 
perature induced by shaking. We have also studied some 
microscopic aspects of density compaction concerning the 
formation of "domains of grains" during the dynamics, 
which we describe in some details. All these results shed 
some light on the origin and the character of the slow 
dynamics in vibrated granular media. An experimental 
check of our findings, which in many cases is still missing, 
would be important to settle our theoretical understand- 
ing of granular media. Actually, the present paper is also 
devoted to compare in some details the two considered 
models and outline their different microscopic properties. 

Before entering the details of the paper, in the follow- 
ing section, for sake of completeness, we briefly explain 
the definition of our lattice gas models and we point out 
their relations to some well known models as the Ising 
model, its Blume-Emery-Griffith counterpart, Edwards- 
Anderson Spin Glasses and driven systems as the Katz- 
Lebowitz-Spohn model. 



II. FRUSTRATED LATTICE GAS MODELS 

The characteristic properties of the Frustrated Lattice 
Gas we introduced to describe the dynamics of dry gran- 
ular materials (in the regime of high packing densities 
and small shaking amplitude), are fully described in Ref.s 
§0. 

The crucial ingredient of these models is, actually, the 
presence of "frustration" in the motion of grains |^,||,|l0| . 
The models consist of a system of elongated particles 
which occupy the sites of, say, a square lattice tilted by 
45''. The particles have an internal degree of freedom 
Si = ±1 corresponding to two possible orientations on 
the square lattice. Nearest neighbour sites can be both 
occupied only if particles do not overlap (i.e., they have 
the right reciprocal orientation) otherwise they have to 
move apart. In our models particles undergo a driven dif- 
fusive Monte Carlo (MC) dynamics. In absence of vibra- 
tions they are subject only to gravity and they can move 
downwards always fulfilling the non overlap condition. 
The presence of vibration is introduced by also allowing 
particles to diffuse upwards with a probability p^p. The 
quantity xq = Pup/Pdown (with pdown = 1 - Pupkas we 
will see, is related to an effective temperature and 
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plays the role of the experimental tap vibrations inten- 
sity. 

Our models can be interestingly mapped on stan- 
dard lattice gas models of Statistical Mechanics, whose 
Hamiltonian has an hard core repulsion term (J — > oo): 
Hhc = JJliij) Ii]{Si, Sj)ninj, where = 0, 1 are occu- 
pancy variables, Si — ±1 are the above spin variables as- 
sociated to the orientations of the particles, J represents 
the infinite repulsion felt by the particles when they over- 
lap. The hard core repulsion function fij{Si, Sj) is or 1 
depending whether the orientations Si,Sj of neighbours 
is allowed or not. 

The choice of fij{Si,Sj) depends on the particular 
model. Here we consider two models: the Tetris and the 
Ising Frustrated Lattice Gas (IFLG). In the Tetris model 
f,,{S,,Sj) is given by: f^f"^{S,,S,) = l/2{S,Sj - 
iij{Si + Sj) + \) here = +1 for bonds along one direc- 
tion of the lattice and = — 1 for bonds on the other. 
This corresponds to a generalised Blume-Emery-Grifhth 
Hamiltonian, and has an "antiferromagnetic" equilibrium 
phase diagram. In order to have a non trivial behaviour 
the dynamics of the Tetris has a crucial purely kinetic 
constraint: particles can flip their "spin" only if many of 
their own neighbours are empty (3, in our simulations), 
as much as in facilitated kinetic Ising models 

Real granular media may have more disorder due to 
broader grain shape distribution or to absence of a reg- 
ular underlying lattice, and each grain moves in the dis- 
ordered environment generated by the others. To de- 
scribe this kind of scenario, we previously introduced a 
model, the Ising Frustrated Lattice Gas (IFLG), made 
of grains moving in a lattice with quenched geometric 
disorder, with the following hard core repulsion function: 
flf'^^iS^.S,) = l/2{e,jS,S, - 1) where = ±1 are 
quenched random interactions associated to the edges of 
the lattice, describing the fact that particles must satisfy 
the geometric constraint of the environment considered 
as "practically" quenched. The IFLG shows a non triv- 
ial dynamics without the necessity to introduce kinetic 
constraints. The Hamiltonian of the IFLG exhibits rich 
behaviours in connection with those of "site frustrated 
percolation" [|7| and Spin Glasses [|6|Jl^. 

The other important contribution to the Hamiltonian 
of a granular media we consider, is gravitational energy: 
Hg = gTlii'^iyi^ where g is the gravity constant and Ui 
is the hight of particle i (we set to unity grains mass 
and lattice spacing). The temperature, T, of the present 
Hamiltonian system (with J = oo) is related to the ratio 
xo = Pup/pdown via the following relation: e^"^^!^ = xq. 
The adimensional quantity F = ln(a:Q ) = T/^, seems 
to play the same role as the amplitude of the vibrations 
in real granular matter 

A further term which could appear in the Hamilto- 
nian is a coupling between neighbouring grains, due, for 
instance, to the presence of some fluids which, due to 
surface tensions, exerts attraction between particles, or 
electrical charges on grains surfaces which might attract 



or repel others grains . Thus the full schematic Hamil- 
tonian may have an other term as the following: 



<ij> 



(1) 



The coupling if 2 must be fixed by the specific interaction 
potential considered. Thus, the complete Hamiltonian 
results to be of the form: 



H — Hhc + He 



H. 



SI 



(2) 



Below, we will generally consider a system of particles 
interacting just via hard core repulsions, i.e., K2 = 0, 
corresponding, for instance, to dry non charged granular 
systems. We will also discuss the presence of attraction 
between neighbouring grains (i.e., K2 > 0), which, we 
show, does not qualitatively change the general scenario 
(at least for K2 not too large with respect to gT). 

It is interesting to notice that the present models are 
very similar to a driven Ising lattice model introduced by 
Katz, Lebowitz and Spohn to describe non equilibrium 
steady states in fast ionic conductors and other systems 
|^,|l|. One of the relevant differences is the boundary 
conditions imposed on particles motions. 
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FIG. 1. The logarithm of the bulk density distribution, 
aP{p), recorded after a random pouring of grains in the sys- 
tem box, as a function of (p — {p)ip)/<J. Here {p)ip is the 
average loose packing density and o its mean square devia- 
tion. MC data from the IFLG (squares) are well fitted by 
a Gaussian, however, data from the Tetris (circles) clearly 
show a violation of such a behaviour. A Gaussian function 
(27r)~^/'^ exp(— a:'^/2) (full line) is shown for comparison. We 
also show an overall Gumbel like fit (dotted line) described in 
the text. 



III. BULK DENSITY DISTRIBUTION 

Monte Carlo simulations of the present models have 
periodic boundary conditions along the horizontal direc- 
tion and rigid walls at bottom and top (here we consider a 
two dimensional sample, but analogous results are found 
in three dimensions as shown in Ref. ) . After fixing the 
tij (which are random in the IFLG), the initial particle 
configuration is prepared by randomly inserting particles 
of given spin into the box from its top and then letting 
them fall down (p„p = 0) until the box is filled. The 
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two basic Monte Carlo moves (the spin flip and particle 
hopping) are done in random order. 

An important quantity to characterise grains packing 
after such a random insertion, is the bulk density distri- 
bution, P{p). Such a quantity is shown in Fig.n^ for the 
IFLG and the Tetris. The data for IFLG concern a lat- 
tice box of size 100 x 100 and for Tetris a box 120 x 1000 
and are averaged over 50000 configurations (the IFLG 
size is smaller since we have also to average over the Cij ) . 

The density where P{p) has its maximum practically 
corresponds to the average loose packing density of the 
system (p)/p {{p)ip ~ 0.739 in the IFLG and {p)ip ~ 0.751 
in the Tetris [Q), whose precise location depends on the 
model and its linear sizes (for instance, boundary effects 
can change it). In Fig.|l| the logarithm of crP{p) is plotted 
as a function of {p — {p)ip)/cr for both data sets [a is the 
density mean square deviation, a ~ 0.002 in the IFLG 
and cr ~ 0.004 in the Tetris). This way to plot P{p) al- 
lows a direct comparison between the two different distri- 
butions. Furthermore, it easily outlines deviations from 
a Gaussian behaviour. Actually, after the above data 
rescaling, a Gaussian variable must have a normalised 
univariate zero mean Gauss distribution function (with 
no adjustable parameters): 



P{x) 



exp(-.'E2/2) 
(27r)i/2 



(3) 



This function is the full line in Fig.Q. Apparently, the 
data for the IFLG model (squares) seem to reasonably 
follow such a Gaussian shape, a fact which well compares 
with experimental data from Ref. p^ . 

The Tetris case, is more controversial from our simu- 
lations. Data from a 100 x 100 sized system seem to fol- 
low the Gaussian distribution found for the IFLG. They 
are not plotted in Figj^ for sake of clarity. In Fig.|l| are 
shown, instead, Tetris data (circles) for the above cited 
largest system size (about 10^ grains) we could simulate 
(their statistics is even better than the one for the smaller 
system). Interestingly, in this case P{p) is not Gaussian. 
For density values below {p)ip^ P{p) rapidly falls off while 
above {p)ip it surprisingly shows an almost exponential 
tail (see Fig.|l|). In Fig.^ we also plot the overall fit ob- 
tained by using a Gumbel-like function: 



P(x) = Po{cxp[— a(6 + x) — e 



(4) 



where Pq ~ 37, a ~ 0.5, b ~ 0.2 and m ~ 4.5. The above 
Gumbel-like distribution is interesting for its relations to 
extreme values statistics but it is difficult to state 
whether such an approximate three parameters fit {Pq is 
fixed by the normalisation) plays an important role in 
the present context. It should be said, however, that ex- 
perimentally a non Gaussian P(£) is usually found after 
long "tapping" of the system [g^l- A finite amount of 
CPU time has not allowed us to go that far. Further 
theoretical and experimental work must be devoted to 
understand the present results. 
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FIG. 2. The presence of attractive interactions among 
grains, may be simulated by turning on K2 in our full Hamil- 
tonian. We show here how such an attraction affects com- 
paction. We plot the bulk density of the Tetris, p, as a func- 
tion of t„, the tap number, for different values of K2 (given in 
units of gV). The system size is 30 x 60 and e~'^^^ = xo = 0.05 
(similar results are found for xo £ [0.01,0.1]). Data are av- 
eraged over 240 taps sequences. The grains attraction makes 
the initial relaxation slower, also if the long times value of 
density is higher the higher is K2- 



IV. COMPACTION IN PRESENCE OF GRAINS 
ATTRACTIVE INTERACTIONS 

After having described the properties of the starting 
random loose packed configuration (for other details see 
Q), in this Section we start discussing the phenomenon 
of compaction. This plays an important role in under- 
standing grains dynamics since it clearly shows the ba- 
sic mechanisms underlying dynamical processes in gently 
shaken granular media subject to gravity 

Referring to recent experiments on compaction dynam- 
ics, here we study the phenomenon of density increase 
due to tapping. In real experiments a "tap" is the shak- 
ing of a container filled with grains by pulsed vibrations 
of given duration and amplitude |^ . In our MC simula- 
tions, in each single tap we apply vibrations of a given 
amplitude a;o to particles for a given, short, duration r 
(i.e., for t e [0,r] we fix x{t) = = const.) Then 
the system is let to find a stationary state for a time 
trepose iu which x{t) = {trepose is choscn to be much 
longer than any relaxation times in absence of shaking 
Q). After each tap we measure the static bulk density 
of the system p{tn) (tn is the n-th tap number). For this 
Monte Carlo experiment with the IFLG model we con- 
sidered a system of size 30 x 60, averaged over 32 different 
eij configurations, and fix r = 32. We refer to Ref. j|] 
for further details on the MC dynamics. 

To describe experimental observations about grain 
density relaxation under a sequence of taps a logarithmic 
law was proposed in Ref. [^: p{tn) = Poo — Apoo/[l + 
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B ln(t„/Ti + 1)]. This law has proved to be satisfied 
very well by relaxation data in the present IFLG model 
as shown in Ref. which can be excellently rescaled 
with experimental data. 

The results from IFLG and Tetris are surprisingly sim- 
ilar ||^,|lo|, but the asymptotic density poo hi the Tetris 
is numerically indistinguishable from 1, thus almost in- 
dependent on xq, a fact in contrast with both IFLG and 
experimental results from Knight et al. . 

In order to have a full comparison of our models with 
a variety of experimental results and to understand the 
effects of small perturbations on the grain interaction 
potential on the general compaction scenario, we also 
dealt with the compaction of systems whose neighbour- 
ing grains feel a finite attraction. This is the case, for in- 
stance, in presence of grains with interstitial fluids, which 
may generally affect the system properties jl^,^ . Thus, 
we turn on K2 in our Hamiltonian eq.(^. 

The results of density compaction of interacting attrac- 
tive grains are shown in Fig.^. The presence of a finite 
K2 doesn't alter the general logarithmically slow features 
of the dynamics in the considered range K2/gT £ [0,8], 
but interesting new phenomena appear. Fig.0 concerns 
the case of the Tetris: the figure shows, for reference, 
the density compaction, p{tn), of a system without inter- 
grains attraction (circles) during a sequence of MC taps 
as those described above; it shows moreover the com- 
paction of other two systems with finite inter-grains at- 
tractions, K2'. the interaction strengths are expressed as 
K2 = 5 (squares) and K2 = 8 (triangles) in units of gT, 
which is the relevant physical energy scale. One can see 
that two seemingly opposite behaviours are present: the 
grains attraction makes the relaxation much slower also 
if the long times value of density is higher the higher is 
K2. For K2 = SfjF, some overall qualitative change in 
the density relaxation curve is also appreciable. 

A. Evolution by "jumps" and "persistence times" 

The density evolution curves shown in Fig.|^ are ob- 
tained as an average over many taps sequences and for 
such a reason they look regular. However, the shape of 
the curve p{tn) in each single run is interesting in itself 
since it is highly irregular and its evolution determined 
by sudden changes. Actually, during several consecutive 
"taps" the systems seems frozen in the same state, then, 
after a new tap, abruptly its bulk density jumps to a new 
value (which only on average is higher that the previous) . 
During such seemingly "dead" time intervals the system 
persists in the same state, with a given bulk density. We 
call "persistence times", tp, the duration of these inter- 
vals. 
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FIG. 3. The logarithm of the "persistence time" distribu- 
tion, P(tp), as a function of tp. tp is the number of elapsed 
taps between two consecutive changes in the bulk density of 
the IFLG during sequences of taps (here with xo ~ 0.0001). 
The lower data (circles) , which are recorded when the system 
density is in the interval [0.735,0.745], and the upper data 
(squares), recorded in the interval [0.745, 0.755], show an ex- 
ponential behaviour (full lines). 

Persistence times have an interesting distribution, 
P{tp), which depends on the system density p. In or- 
der to increase the accuracy, data plotted in FigJ| are 
values obtained after averaging in a given density inter- 
val. At low densities (close to {p)ip), P{tp) certainly has 
exponential shape: 

Pitp) ^ e-*-/*p (5) 

In the interval, 0.735 < p < 0.745, we find t° ~ 0.40. 
In a higher densities interval 0.745 < p < 0.755, P{tp) 
seems to becomes broader (see Fig.||). It is extremely 
difficult to collect sufficiently clean data for high values 
of tp in the high density region (due to the required very 
long CPU times) in order to investigate such a "broad- 
ening " of the -P(0- An exponential fit (with t° ~ 0.42) 
is still possible [g^. These results are plotted in Fig||, 
where we show data for the IFLG (size 30 x 60) recorded 
during sequences of taps with xq — 0.0001 {K2 — 0). 
Similar behaviours of P{tp) are observed in the region 
Xq G [10~^, 10""'^], showing that not only the overall prop- 
erties of the dynamics do not change in such a "shaking" 
interval. 

We discuss below, from a microscopical point of view, 
how density changes are made possible by the local reor- 
ganisation of the structure of grains clusters. 
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FIG. 4. In our lattice systems, we define a "cluster" as 
a set of nearest neighbour grains. In the present figure (con- 
cerning a lattice of size 30 x 60), we show the number of 
clusters, Ndust, as a function of the taps index i„ during a 
MC sequence of taps in both Tetris (filled circles) and IFLG 
(empty squares). In both kind of models are there very few 
clusters for unit length and their number very slowly decreases 
with t„. The data are averaged over 100 taps sequences with 
xo = 0.1. 



V. DOMAINS OF GRAINS IN THE PACK 



^4 

(0 

I 

E 

V) 






OTetris 




□ IFLG 



o o 



10^ 



10 



10' 



10' 



0.015 






OTetris 


□■■ 


□ IFLG 



0.010 



a 
a 



0.005 - 

D 



Exploiting the microscopic character of our models, we 
have access to the details of the above compaction process 
and we explore the grains packing structure. 

In our lattice systems we can define a "cluster" of 
grains as a set of nearest neighbouring particles. Due 
to grains shapes incompatibilities, empty sites can be 
present in the lattice pack. For instance in the Tetris, 
clusters can be pictorially described to be mainly thick 
"vertical" clumps separated by almost linear sequences of 
holes. Actually, the above definition of "clusters" must 
be refined to be applicable to a real granular medium, 
but on a lattice it is the most natural. 

In the randomly prepared initial state are there typi- 
cally very few clusters per unit length (in a 30 x 60 sized 
lattice one has a total of 3-4 clusters, see Fig.^. By tap- 
ping, the number of clusters, Ndust, approximately loga- 
rithmically decreases with the taps index tn (see Fig.^. 
In the same time the larger cluster grows and the others 
shrink (see Fig.^). The data shown in Fig.^ and Fig.|^ 
refer to tapping sequences with "amplitude" xq = 0.1, 
and analogous results were found for xq = 0.01. 
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FIG. 5. Top: we plot the difference of the average sizes of 
the largest cluster in the system, Sm, minus the same quan- 
tity after the first tap, 5'A/(ti), as a function of the taps index 
t„, during a MC taps sequence with 2:0 = 0.1 in both Tetris 
(circles) and IFLG (squares). The box size is 30 x 60. The val- 
ues at ti are: SM{ti) = 919 in the IFLG, and SM{ti) = 687 
in the Tetris. Bottom: for a clear comparison, we also plot the 
bulk density variation of the same systems during the same 
taps sequences. The values at ti axe: p{t\) = 0.743 in the 
IFLG, and p{t\) = 0.830 in the Tetris. Data are averaged 
over 100 taps sequences. 

Also if the overall behaviour is similar, the compaction process 
of the two considered models is microscopically very differ- 
ent. In the Tetris it is originated by the growth of the largest 
cluster at expense of the smaller ones, as in spinodal decom- 
position. In the IFLG on the contrary, the growth of the 
largest cluster is extremely weak: the compaction is mainly 
due to the "expulsion" of holes from the largest clusters and 
the consequent optimisation of grains arrangement. 



The analysis of cluster properties reveals interesting 
facts about the compaction mechanisms and outlines a 
main difference between the compaction in IFLG and 
Tetris: it concerns the relative size of the few present 
clusters. In the IFLG the second largest cluster is, since 
the beginning, made of very few grains (0(1)); in the 
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Tetris, the first and second largest clusters are instead, 
at the beginning, almost of the same order of magnitude. 

During the taps sequence, in the Tetris the largest clus- 
ter grows at expense of the smaller ones. As shown in 
Fig.H, after 10'^ taps, a 1% increase of the bulk density, 
p{tn), is in direct correspondence with a 1% increase 
in the largest cluster mass, SM{tn)- This mechanism 
strictly recalls ordered domains growth in spinodal de- 
composition As a matter of fact, the Tetris is a 
"kinetically constrained" antiferromagnet and our clus- 
ters explicitly correspond to the definition of the "Fisher 
droplets" in the system In the IFLG, on the con- 
trary the growth of the largest cluster is extremely weak 
(see Fig.||) since the other few clusters size is already of 
the order of very few grains: now a 2% increase in the 
bulk density corresponds to only a 0.1% increase in the 
largest cluster mass. In the latter case the compaction 
of the system is, thus, mainly due to the "expulsion" of 
holes from the largest cluster and the consequent optimi- 
sation of grains spatial arrangement. 

In the present perspective the basic mechanisms of do- 
mains growth under gravity in the IFLG and Tetris seem 
to be very different. To understand what scenario is 
closer to reality, if any of the present, an experimental 
investigation of this kind of properties would be very im- 
portant. 

The description of compaction with a diffusion equa- 
tion with a density dependent non linear diffusion coeffi- 
cient was proposed in Ref.s It resembles analogous 
behaviours found close to the glass transition in glass for- 
mers [^. 



granular packs During such a process, the dynam- 
ics becomes slower and slower since grains self-diffusivity 
decreases with density Q. At high densities the sys- 
tem approaches a structural arrest at a "jamming" tran- 
sition point, where self-diffusivity becomes zero, a fact 
very similar to the "freezing transition" in glass formers 
pT|-p^. Interestingly in the IFLG, the structural arrest 
due to self-diffusivity suppression coincides with a spin 
glass transition Q . 

The microscopic nature of the models allows the study 
of the properties of their internal packing structure, as 
grains clusters, which give a deeper understanding of the 
basic mechanisms underlying compaction. The study of 
"grains clusters" shows that in the IFLG model com- 
paction originates from the reduction of the number of 
holes in the only present large cluster by optimisation of 
grains spatial arrangement; in the Tetris, instead, com- 
paction stems from the growth of the largest cluster mass 
itself at expense of the smaller ones, as in spinodal de- 
composition. Actually, this result points towards the ne- 
cessity of a better understanding of the real long time 
dependence of compaction in the Tetris and IFLG. 

Many of the results presented here are still waiting for 
an important experimental verification. 

Work partially supported by the TMR Network Con- 
tract ERBFMRXCT980183 and MURST-PRIN 97. 



VI. CONCLUSIONS 

Summarising, the central body of the present paper 
has dealt with the dynamical behaviours of two frustrated 
lattice gas models (the IFLG and Tetris) introduced to 
describe gently vibrated granular media HJl^. These 
microscopic models are characterised by a gravity driven 
diffusive dynamics where the basic ingredient is the pres- 
ence of geometric frustration in particles rearrangement. 
Interestingly, although granular media are non-thermal 
dissipative systems, the present models can be casted in 
the Hamiltonian formalism of standard lattice gases of 
Statistical Mechanics in presence of an effective temper- 
ature (different from zero during shaking). 

To characterise the grains packing structure we have 
recorded the loose packing density distribution function, 
P{p), which in the case of the IFLG well compares with 
experimental results. More interesting, but also more 
controversial is the non Gaussian P{p) recorded in the 
Tetris. 

Also in presence of attractive surface grains interac- 
tions, the models exhibit a slow compaction when subject 
to gentle shaking in presence of gravity, a compaction 
extremely close to what is experimentally observed in 
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